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The BES collaboration has recently observed a strong enhancement close to the proton- 
antiproton, pp threshold in the J/if> decays into -ypp. Such a structure can be explained 
by a traditional nucleon-antinucleon, NN, model. The near threshold 11 5b bound state 
and/or the well-established 13 Po resonant state found in this NN interaction can ade- 
quately describe the BES data. 
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1. Introduction 

Existence of near threshold bound states or resonances in nucleon-antinucleon, NN, 
interaction is a challenging matteffl Low-energy scattering could indicate the pres- 
ence of such structures by determining the scattering lengths for 2I + 1 - 2S + 1 t,j states. 
Here I denotes the isospin (0 or 1), S the spin (0 or 1), L the angular momentum 
and J the total angular momentum. An alternative is to use formation experiments. 
At the Beijing electron-positron collider, the BES collaboration has observed a 
resonant-like behavior in the pp invariant mass spectrum from radiative J/ip — > "fPP 
decays^. The present work studies the physics of slow pp pairs produced in J/ip 
decays, using J PC conservation, P being the parity and C the charge conjugation. 
Here we rely on the Paris NN potential model. 



2. Close to Threshold Proton- Antiproton Final State Model 
2.1. The low-energy nucleon-antinucleon interaction 

The Paris NN interaction is built up from a state dependent optical potential. 
The long range, r > 1 fm, real part is obtained by G-parity transformation of the 
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Paris NN potential, the two-pion exchange of which is calculated via dispersion 
relations from pion-nucleon scattering data. The short ranges, r < 1 fm, real part 
and absorptive part, with a form suggested by calculation of NN annihilation into 
two mesons or resonances, are both determined through fit to the NN data. In 
the different versions, the short range parameters are readjusted by fitting to new 
data. The Paris 82 potential, fitted to pre-LEAR (CERN) data, mainly elastic 
pp (isospin 1 + isospin 0) data , has a x 2 /data of 2.8 for 915 data. The Paris 94 
potential uses LEAR data, in particular pp — > nn (7 = 1 — 1 = 0) data, and has 
a x 2 /data of 2.46 for 3295 data. In the Paris 99 versio more recent LEAR data, 
in particular for pp — > fin, were used leading to a x 2 /data of 2.95 for 3814 data. 
The Paris 04 models is constrained by fitting to the 1999 data plus the scattering 
lengths extracted from antiprotonic hydrogen and deuterium dat£0 and to the total 
np cross-section^. It has x 2 /data = 3.19 for 3934 data. 

2.2. Allowed slow pp final states 

The J PC conservation (J PC = 1~~ for J/ip) limits the number of slow pp final 
states. These correspond to pairs of small M pp — 2m p with M pp being the invariant 
pp mass and m p t he p roton mass. The allowed states are listed in Table 1. Some two- 
particle analogues^ are listed in the second column. The last column indicates the 
relative angular momentum between 7 or ir and the pp pair h. The BES experiment 
angular distribution favors a pseudoscalar 1 5'o or a scalar 3 Pq h final state. 



Table 1. The slow pp pairs states permitted in the radiative J /if) — > 7pp. 



decay mode 


analogue 


J PC [ 7 or tt] 


J PC [pp] 


h(pp) 


relative i 


7PP( 1 S'o) 


jr) (1444) 


1 — 


0-+ 


pseudoscalar 


1 


7Pp( 3 P ) 


Y/o(mo) 


1 — 


0++ 


scalar 





7Pp( 3 Pi) 


7/1 (1285) 


1 — 


1++ 


pseudovector 






2.3. Specific final-state interaction model 

The transition amplitude from a channel i to a channel /, in a multichannel system 
at low energy described by a S-wave K matrix, can be written as Tjf — Aif(l + 
iqfAf . Here Aif is a transition length, Aff the scattering length in the channel 
/ and qf the momentum in this channel™ The / channel scattering amplitude can 
also be expressed as Tjf — Aff(l + iqjAff)^ 1 . For a P wave close to threshold, 
Aff = A?f Of and Aif = Af^qf where A P j is the scattering volume. Up to terms in 
q) one has T if = (A lf /A f f)Tff = CT ff /q p = Ct L . The quantity C = A if qf/A ff 
represents the unknown formation amplitude and |£l| 2 = \Tff/q p | 2 is the final state 
interaction factor in a given pp partial wave. In terms of the phase shifts 8l and 
inelasticities 7]l of a given NN interaction one has — (j/ie 2 " 5 ' — l) /(2iqj L+1 ). 
The function C is parametrized by |C(a;)| 2 = qf(co + c\x) where x = M pp — 2m p 
and qf = [x{m p + cc/4)] 1 / 2 . 
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3. Results and Conclusions 

The final state interaction factors |irj 2 for the 1 Sq and 3 Pq states and for the 
different versions of the Paris NN are compared to the BES data-in Figs, f and 
2. The Co and c\ parameters are determined by requiring |T"i/| 2 of Paris 04 to be 
close to the events distribution as given in Fig. 3 of Ref. 2 at x — 7 MeV and 
x = 66.2 MeV. For ^o, c = 1.18599, c x = 0.00299 and for 3 P , c = 2.5206, 
ci = 0.0269. As seen in Fig. 1, the data is well reproduced by the Paris 04 NN 
interaction. This interaction has a 11 S'o bound state located at x ~ —4.8 MeV and 
with a width T of 52.5 MeV. Paris 99 has also a bound state at x = —69 MeV with 
r = 46 MeV. There are no bound states for Paris 94 or Paris 82. All Paris models 
have a 13 Po resonance of mass ~ 1876 MeV and r ~ 10 MeV. They all reproduce 
the near threshold BES enhancement as seen in Fig. 2. 

In conclusion, the near threshold pp enhancement seen in BES collaboration^ 
can find a natural explanation from a traditional model of pp interaction. The 11 So 
bound statJ^ needs confirmation. The well established 13 Pq resonance originates 
from the strong attraction of the one-pion exchang Each of these states gives 
a reasonable representation of the BES radiative J/ip — > "/pp decay data. They 
correspond to the S or P wave Breit Wigner resonance functions considered by 
BES collaboration in their fit to the datcPl 



J/V^ypi 



Paris 04 (Kef. 6) 

Paris 99 (Ref. 5) 

Paris 94 ( Rcf 4 ) 

Paris Ha (Kef. 3) 

|t(x)| 2 =|T(x)| 3 /(q(c o +c iX )) 
x=M _-2m , q=(x(m +x/4)) 1/s 
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Fig. 1. The 1 So final state factor compared to BES datsPl 
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